Orientational order and glassy states in networks of semiflexible polymers 
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Motivated by the structure of networks of cross-linked cytoskeletal biopolymers, we study orienta- 
tionally ordered phases in two-dimensional networks of randomly cross-linked semiflexible polymers. 
We consider permanent cross-links which prescribe a finite angle and treat them as quenched disor- 
der in a semi-microscopic replica field theory. Starting from a fluid of un-cross-linked polymers and 
small polymer clusters (sol) and increasing the cross-link density, a continuous gelation transition 
occurs. In the resulting gel, the semiflexible chains either display long range orientational order or 
are frozen in random directions depending on the value of the crossing angle, the crosslink concen- 
tration, and the stiffness of the polymers. A crossing angle 6 ~ 2n/M leads to long range M-fold 
orientational order, e.g., "hexatic" or "tetratic" for 6 = 60° or 90°, respectively. The transition to 
the orientationally ordered state is discontinuous and the critical cross-link density, which is higher 
than that of the gelation transition, depends on the bending stiffness of the polymers and the cross- 
link angle: the higher the stiffness and the lower M, the lower the critical number of cross-links. 
In between the sol and the long range ordered state, we always expect a gel which is a statistically 
isotropic amorphous solid (SIAS) with random positional and random orientational localization of 
the participating polymers. 

PACS numbers: 87.16Ka,82.70Gg,61.43Er 
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I. INTRODUCTION 



The cytoskeleton is a network of linked protein fibers 
which plays an important role in several functions of eu- 
karyotic cells such as maintenance of morphology, me- 
chanics and intracellular transport [1]. Cytoskeletal 
fibers, such as F-actin can be described as semiflexible 
polymers with a behavior intermediate between the two 
extreme cases of rigid rods and random coils. The func- 
tion of the actin cytoskeleton is modulated by a large 
number of actin-binding proteins (ABPs) [2, 3]. The or- 
ganization of actin filaments into networks is regulated by 
ABPs which can be classified into two broad categories; 
cross-linkers promote binding of the filaments at finite 
crossing angles whereas bundlers promote formation of 
bundles consisting of parallel or antiparallel filaments. 
F-actin is a polar semiflexible polymer and some ABPs 
bind filaments with a specific polarity whereas others are 
not affected by the filament's polarity. In order to un- 
ravel the physics of these complex aggregates, in vitro 
solutions of actin filaments with controlled cross-linkers 
have been studied [4]. 

The stiffness of semiflexible filaments gives rise to 
orientational correlations and promotes the formation 
of structures with long range orientational order. The 
isotropic-nematic transition in solutions of partially flex- 
ible macromolecules has been studied theoretically us- 
ing the Onsager approach [5] or the Maier-Saupe ap- 
proach [6, 7]. In Ref. [7], the role of the solvent is 
taken into account and a very rich phase transition ki- 
netics is predicted. Aggregation and orientational or- 



dering of Lcnnard- Jones macromolecules with bending 
and torsional rigidity, with or without solvent, has re- 
cently been investigated in Ref. [8] using molecular dy- 
namics simulations. Experimental investigations of the 
isotropic-nematic transition in lyotropic F-actin solutions 
have been carried out in Ref. [9] and measurements of 
the associated order parameter are presented in [10]. 

The excluded volume effect is not the only mechanism 
which drives an isotropic-nematic transition in solutions 
of semiflexible polymers. Assemblies of cytoskeletal fil- 
aments exhibit a structural polymorphism due to inter- 
actions mediated by a wide range of ABPs, as shown in 
Ref. [11]. Theoretical attempts to study the formation of 
ordered structures in this kind of systems involve a gener- 
alized Onsager approach [12-14], a Flory-Huggins theory 
[15], and a semi-microscopic replica field theory [16]. In 
Rcfs. [12-15] the filaments are modeled as perfectly rigid 
rods whereas in Ref. [16] they are considered to be semi- 
flexible polymers and the thermal bending fluctuations 
(finite persistence length) are fully taken into account. 

In vitro structural studies by Wong et al. [17] of F- 
actin in the presence of counterions have shown the for- 
mation of sheets ( "rafts" ) with tetratic order without any 
cross-linking proteins. It appears that electrostatic in- 
teractions are the main mechanism behind the effective 
"7r/2" cross-linking of the actin filaments in this experi- 
ment. A theory for tetratic raft formation by rigid rods 
and reversible sliding ll ir/2" cross-linkers has been pro- 
posed by Borukhov and Bruinsma [13]. 

In Ref. [16], a three-dimensional melt of identical, 
fixed-contour-lcngth semiflexible chains is considered and 
random permanent cross-links are introduced. The cross- 
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links are such that they fix the relative positions of the 
corresponding polymer segments and constrain their ori- 
entations to be parallel or antiparallcl to each other. The 
aim of the present study is to examine the effect of cross- 
links which prescribe a finite crossing angle on a two- 
dimensional version of the previous model. We distin- 
guish the case of sensitive cross-links that are perceptive 
of the polymers' polarity and unsensitive cross-links that 
are not. In contrast to previous studies where the rigid 
rods have a finite width which allows for nematic order- 
ing a la Onsager, our semiflexible polymers are considered 
one-dimensional objects and the sole cause for the emer- 
gence of oricntationally ordered phases is the interplay 
of the finite persistence length of the polymers and the 
cross-link geometry. 

Using the polymer stiffness and the cross-link density 
as control parameters, we obtain a phase diagram which 
involves a sol and various types of oricntationally ordered 
gels. For appropriate values of the control parameters 
and crossing angle 8 ~ 2n/M;M £ Z, we predict the 
emergence of an exotic gel with random positional lo- 
calization and M-fold orientational order (e.g., hexatic 
for M — 6 or tetratic for M — 4). Similar phases have 
been predicted [14] for a different system in three space 
dimensions: a semi-dilute solution of charged rods with 
finite diameter in the presence of polyvalent ions that 
have the function of (non-permanent) cross-linkers and 
may favor various crossing angles. Besides the orien- 
tationally ordered phase, we also predict a statistically 
isotropic amorphous solid (SIAS) with random positional 
and orientational localization of its constituent polymers 
appearing right at the gelation transition. 

This article is organized as follows. In Section II, we 
introduce our model and the Deam-Edwards distribu- 
tion which parametrizes the quenched disorder associated 
with the cross-links. In Section III, the disorder-averaged 
free energy is presented as a functional of a coarse- 
grained field which plays the role of an order parame- 
ter. Then, in Section IV, we discuss different variational 
Ansatze that express positional and orientational local- 
ization. The symmetries imposed by the cross-linking 
constraints allow the emergence of specific oricntation- 
ally ordered phases. The corresponding free energies are 
calculated variationally in the saddle-point approxima- 
tion and a phase diagram is obtained. We summarize 
and give an outlook in Section V. Finally, details of the 
calculations are given in the appendices. 



II. MODEL 

We consider a large rectangular two-dimensional vol- 
ume (area) V which contains N identical semiflexible 
polymers. A single polymer of total contour length L 
is represented by a curve in 2d space with r(s) denoting 
the position vector at arc-length s. Bending a polymer 
costs energy according to the effective free energy func- 
tional ( "Hamiltonian" ) for the wormlike chain (WLC) 



model [18] 

n WLC (Ms)}) = l^^as (^) ■ (1) 

Here we have introduced the tangent vector t(s) = 
dr(s)/ds and chosen a paramctrization of the curve, in 
accord with the local inextensibility constraint of the 
WLC, such that |t(s)| = 1. The position vector is re- 
covered from t(s) as r(s) = r(0) + / t(r)dr. Hence 
the conformations of a single polymer, that can be bent 
but not be stretched, are alternatively characterized by 
t( T )> < r < L, and r(0). The bending stiffness is 
denoted by k; it determines the persistence length L p 
according to k = L p kBT/2. Throughout the rest of the 
paper we set UbT — 1. The WLC model can describe stiff 
rods, obtained in the limit L/L p — > and fully flexible 
coils, obtained in the limit L p / L — > 0. 

The mutual repulsion of all monomers is modeled by 
the excluded- volume interaction 

N 2 

^ ev = W X/ X/ ^f k l.wl/ 5 k,m| 2 (2) 
k^O rogZ 

with 

1 N L 

Pk, m := E T / ds c?kr,(s) &imM8) ( 3 ) 

i=l ^° 

being the Fourier transformation of the positional- 
orientational density. tj(s) = (cos ^(s), sin ^(s)) de- 
notes the orientation of monomer s on polymer i. The 
coefficients A 2 k | m depend only on the absolute value of 
the vector k in order to preserve the rotational symme- 
try of the system. They are later chosen large enough 
to provide stability with respect to density modulations. 
For details, see Appendix C. 

The Hamiltonian T-Lwlc is invariant with respect to in- 
terchanging head and tail of the filaments, i.e., the energy 
is unchanged under reparametrizations of the contour of 
one polymer i by {r,(s) — > Yi(L — s),Vs £ [0,L]}. How- 
ever, in the following, we want to consider filaments with 
a definite polarity which could for example arise due to 
the hclicity of F-actin. The WLC Hamiltonian is not sen- 
sitive to such a polarity, but the cross-links may or may 
not differentiate between the two states of the filament, 
as discussed below. 

We now introduce M permanent (chemical) cross-links 
between pairs of randomly chosen monomers. A single 
cross-link, say between segments s and s' on two chains 
i and j, constrains the two polymer segments to be at 
the same position, i.e. r^(s) = rj(s'), and fixes their rela- 
tive orientation by a constraint expressed by the function 
A (tj(s), tj ;(s'), 6). The constraint of relative orientation 
is most easily formulated in polar coordinates where A 
just fixes the crossing angle tp — ip' to some prescribed 
value 0. In the case of cross-links sensitive to polarity we 
simply have 

A^-^',0) = Sty-tf-O), (4) 
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where the weight of the delta function is such that 

/o" Itt ^WO = ^ n ^ ne unsensitive case on the other 
hand, changing head and tail of either filament leads to 
four equivalent situations that correspond to two differ- 
ent crossing angles (see Fig. 1). In this case the cross-link 
(4) thus appears with the two equally likely cross-linking 
angles 9 and 9 + it. We point out that in our model the 
polarity is only recognized by the cross-links and does 
not enter into the Hamiltonian. 



Here the thermal average < ... > n is taken with respect 
to the Hamiltonian % := Hwlc +Hev and the partition 
function depends on the cross-link realization C. The 
free energy F is expected to be self-averaging so that we 
are allowed to compute the disorder averaged free energy 
[F] = — [InZ], where [...] denotes the "quenched" average 
over all cross-link configurations according to some dis- 
tribution V{C). We assume that the different realizations 
obey a Deam-Edwards-likc [19] distribution 




FIG. 1: 4 possibilities for two filaments to be joined by an 
unsensitive cross-link. 



The delta function in (4) is the simplest way to model 
the angular constraint of the cross-links. It is however 
much more realistic to introduce an effective angular 
cross-linking potential that allows for fluctuations around 
the preferred direction. A simple model for these "soft" 
cross-links is given by 



3 7 cos(t/i— V>'- 



(5) 



The cross-link stiffness parameter 7 controls the variance 
of the fluctuations of the angle and it is clear that in the 
limit 7 — > 00 we recover the simple delta function model 
of "hard" cross-links. In the following we will first explore 
the case of hard cross-links, because it is mathematically 
simpler. It gives rise, however, to artifacts which dis- 
appear when considering the more physical model which 
favors certain angles but does not enforce them strictly. 

The cross-links are permanent and do not break up or 
rebuild. Hence we are led to study the equilibration of 
the thermal degrees of freedom {r, (s)} in the presence of 
quenched disorder represented by a given cross-link con- 
figuration C = {i e , j e ; s e , s' e }e£i which is characterized 
by the set of pairs of polymer segments that are involved 
in a cross-link. The central quantity of interest is the 
canonical partition function 
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Z(C) = ( 1] S (*M ~ r je (s' e ))A (t ie (Se) - tM),e) 



M 



M 



H 



\8(vM-v je {s' e )) 



(0) 



v{Cm) k m {— 

\e=l / 

(7) 

implying that polymer segments that are likely to be close 
to each other in the uncross-linked melt also have a high 
probability of getting cross-linked. Note that the proba- 
bility of cross-linking is independent of the relative orien- 
tation of the two segments, so that the cross-link actually 
reorients the two participating segments. The parameter 
/x 2 controls the mean number of cross-links in the system: 
[M]/N ~ 0( M 2 ) . 

Given the Hamiltonian, the constraints due to cross- 
linking and the distribution of cross-link realizations, the 
specification of the model is complete and we proceed to 
calculate the disorder averaged free energy [F]. 



III. REPLICA FREE ENERGY 



The standard method to treat the quenched disorder 
average is the replica method, representing the disorder 
averaged free energy as [F] = —[In Z] = — lim„_j.o([^"] — 
l)/n, i.e., we have to calculate the simultaneous disor- 
der average of the partition sums of n non-interacting 
copies of our system. In the end, we extract the disorder 
averaged free energy [F] from [Z n ] as the linear order 
coefficient of its expansion in the replica number n. 

In this section, we present only the essential formulas 
and present the details of the calculation in Appendix B. 
The disorder average, [Z n ], gives rise to an effectively uni- 
form theory, however, with a coupling of different repli- 
cas: 



\Z n \ cx ( exp 



fi 2 V 
2N 



N 



]T / <KM*)-W))A(&(*)-&(*V) 



H 



n+1 



(8) 



where the average is over the (n + l)-fold repli- catcd Hamiltonian T~L. To simplify the notation we 
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have introduced the abbreviations f = (r°, r 1 , . . . , r n ) , 
tfj = (tp 1 , . . . , tp n ) , and the shorthand notation J = 

(l/L)/ L ds. Furthermore, 5{r) = nLo^( r 1 and A 
denotes for sensitive cross-links 

n 

A,(&0) = Y[S(^-9) (9) 

Q = l 

and 

if" " 1 

A u $, o) = - n <5(^ Q - e) + n w - ( + ^)) \ 

[a=l a=l J 

(10) 

if they are unsensitive to polarity of the filaments. For 
soft cross-links we introduce the corresponding defini- 
tions. 

Different polymers are decoupled via a Hubbard- 
Stratonovich transformation, introducing collective 
fields, fi(x, <p), whose expectation values are given by 

N 

<fi(x^)> = -£ / <*(*-*i(*)W0 -&(*))>• (ii) 
i=i ■ Js 

The field (Q) quantifies the probability to find monomer 



s on chain i at position x° in replica 0, at position x 1 in 
replica 1,... and at position x ra in replica n and similarly 
to find it oriented in the direction e 1 = (cos(< / 9 1 ), sin((/? 1 )) 
in replica 1,... and oriented in the direction e" = 
(cos(</? n ), sm((p n )) in replica n. Sometimes it is conve- 
nient to use its equivalent representation in Fourier space 

i N r - 

(ft(k,m)) = -£ / (c^M e***M>. (12) 
i=i • /s 

In terms of these collective fields the effective replica 
theory is given by 

[z n ] ~ Jv{n}^ v {-NTm))). (13) 

The replica free energy per polymer reads 

:F =4^^ m2 - inzmh (i4) 

m 

with the effective single chain partition function 



*({0}) 



cxp f ^-Jlk 51 A n / cxp { ~ ^ f ( s ) } cxp { ~ i™^ 3 ) } 



V21/ 



(15) 



The average (...)^5[ c refers to the (n + l)-fold replicated 
Hamiltonian of a single wormlike chain. The sum yV 
over replicated wave vectors is restricted to the com- 
bination of the so-called "0-replica sector" (ORS) that 
contains only the point k = (0, . . . , 0) and the "higher- 
replica sector" (HRS) where at least two wave vectors in 
different replicas are non-zero: k a ^ and k^ ^ with 



Note that in the above overview we left out the "1- 
replica sector" (IRS) that consists of (n + l)-fold repli- 
cated vectors k Q where only one entry (the ath) is non- 
zero, i.e. k Q = (0, k", 0), and rh being arbi- 
trary. The corresponding fields 51 Q (k, rh) describe regular 
macroscopic density fluctuations (modulated states). In 
this work we focus on the properties of the macroscopi- 
cally translationally invariant amorphous solid state and 
assume that the inter-polymer interactions (2) are such 
that periodic density fluctuations are suppressed. Sec 
Appendices B and C for more details on that issue. 



IV. VARIATIONAL APPROACH 

We shall only discuss the saddle-point approximation 
to the field theory on the right hand side of (13) replacing 
the field 51 by its saddle point tt sp which has to be calcu- 
lated from 6J-/5Cl\n s = 0. In fact, even the saddle point 
equation is too hard to solve, because the conformational 
distribution of the WLC is very complex. Similarly, it is 
not possible to perform a complete stability analysis of 
the Gaussian theory. Hence we restrict ourselves to a 
variational approach and in the following we are going 
to construct Ansiitze which capture the symmetry of the 
different physical states which may emerge. 

What behavior do we expect? Upon increasing the 
number of cross-links up to about one per polymer, there 
should be a gelation transition from a liquid to an amor- 
phous solid state. A finite fraction of the polymers form 
the percolating cluster and are localized at random po- 
sitions. The other polymers belong to finite clusters or 
remain un-cross-linked and are still free to move around 
in the volume. This scenario has been found to be valid 
for a variety of different models [20-24] . A similar replica 
field-theoretic approach has beeen used by Panyukov and 



FIG. 2: Sketch of the tetratic phase (M=4). 



Rabin to study the properties of well cross-linked macro- 
molecular networks [25, 26]. 

For scmiflcxiblc polymers the positional localization in 
the macroscopic cluster is accompanied by orientational 
localization. The cross-links create locally an orienta- 
tional structure. Supposing that the semiflexible poly- 
mers are rather stiff, a long range orientationally ordered 
state can be established. Otherwise we may find an ori- 
entational glass, i.e., the directions of the polymer seg- 
ments are frozen in random directions in analogy to the 
low temperature phase of a spin glass[27]. We call such 
a phase statistically isotropic amorphous solid (SIAS). 

Can we expect long range orientational order for ar- 
bitrary crossing angles 91 We first consider a special 
case, namely that the crossing angle is such that an 
integer multiple of the crossing angle, M9, equals 27r, 
e.g. 9 = 120° and M = 3. Such a choice of crossing 
angle allows for orientationally ordered states provided 
the chains are sufficiently stiff. In the case of sensitive 
cross-links we expect M-fold discrete rotational symme- 
try, in the case of unsensitive cross-links and odd M we 
expect 2M-fold symmetry because cross-links are estab- 
lished including both the angles 9 and 9 + ir (see Fig. 1). 
Sketches of gels with long range four- fold order (9 = 90°) 
or long range three-fold or six-fold order respectively 
(9 = 120°/ 60°) arc shown in Fig. 2 and Fig. 3. The 
case 9 = kjj with k = 1 , . . . , M — 1 leads in a similar 
fashion to an M-fold symmetric phase. 

Let us now consider the constraints on the order pa- 
rameter imposed by symmetries. To this end we note 
that the Hamiltonian (1), the cross-link constraints (6) 
and the Deam-Edwards distribution (7) arc invariant 
with respect to uniform translations and rotations of all 
particles positions {rj(s)}. The disorder averaged parti- 
tion function of Eq.(8) is invariant under spatially uni- 
form translations and spatially uniform rotations of each 
replica separately. 

Only the fluid state is expected to have the full sym- 
metry of the partition function Eq.(8) because here, the 
polymers and finite clusters of polymers are free to sam- 
ple the complete volume and can take any orientation. 



FIG. 3: Sketch of the triangular phase (M=3) for sensitive or 
the hexatic phase (M=6) for unsensitive cross-links. 



This implies for the order parameter 

(a(k,rh)) = 6 ifi S^ d , (16) 

i.e., each replica is separately invariant under translations 
and rotations. 

We will only investigate gels with random localization 
of a fraction of the particles. In other words, we re- 
strict ourselves to amorphous solids and do not consider 
periodic density fluctuations. It might be of interest to 
study one-dimensional periodic density modulations with 
an overall orientation of the WLCs perpendicular to the 
wave vector,- reminiscent of bundles. However this is 
not the topic of the present paper, where we stick to the 
incompressible limit. 

If the translational and rotational symmetry is broken 
spontaneously, then we expect a non-trivial expectation 
value of the local density 

p, )S (x, ip) := 5(x - Ti(s)) 5(<p - i>i(s)) (17) 

in a particular equilibrium state. For example, in the 
amorphous solid phase a finite fraction of particles should 
be spatially localized at preferred positions and possi- 
bly oriented along preferred directions. A simple Ansatz 
quantifying such a scenario is the following: 

( Pi , s (x, p)) OC c-^-^Vos^-cM ( 18) 

The vector a is the preferred mean position of monomer 
s and the thermal fluctuations around the preferred 
position are controlled by the localization length £. The 
preferred orientation is given by (po and r\ parametrizes 
the variance of the orientational distribution. 

An amorphous solid is macro scopically translational 
invariant, so that the spontaneous symmetry breaking 
of translational invariance is local. All macroscopic ob- 
servables should display translational symmetry and in 
particular all moments of the local density 

n J (ki,...k J> m 1 ,...ro J ) (19) 

: = ^E^IiX 5 (^, s ( k i 5 ™i)>-(^, s ( k J>™j)) 

~ <5ki + ...k 7 ,0 
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are non-zero only, if the wave vectors add up to zero (see 
reference [21] for a more detailed discussion). 

As far as the rotational symmetry is concerned, we will 
consider two possibilities: a statistically isotropic state 
(SIAS) as well as states with true long range orienta- 
tional order. In the former case the spontaneous sym- 
metry breaking of rotational invariance is local and re- 
stored globally analogously to what happens in the case 
of translational symmetry. All macroscopic properties, 
such as the moments, are non-zero only, if the angular 
momentum "quantum" numbers add up to zero: 

fi J (k 1 , ...kj, mi, ...mj) ~ <5ki+...kj,o <W+...mj,o ( 20 ) 

For the long range ordered case the simplest Ansatz 
consists in assuming a local M-fold symmetry. This im- 
plies for the moments 

J 

n J (ki,. ..kj, mi,. .. mj) ~ <Sk 1 +...k. /l o JJ <W,zm , (21) 

a=l 

i.e. each m a has to equal an integer multiple of M. 

Note that a rotation affects both rj(s) and t,(s) and 
that consequently, a system with e.g. M-fold symmetry 
should be described by an Ansatz Q(k, ip) which is 
only invariant under common rotations of the k" and 
ip a . We have chosen our approach for simplicity and 
incorporated the symmetry with respect to the orien- 
tational and position vector separately. Consequently, 
the Ansatz is symmetric under individual rotations 
of the spatial and angular argument. An interesting 
extension of the present work would involve the study 
of Ansatze which couple positional and orientational 
degrees of freedom. Some preliminary work in this 
direction has been done in the context of random net- 
works of covalently connected atoms or molecules [28, 29]. 



Rephrasing our results in replica language, frozen 
fluctuations in a single equilibrium state correspond to 
a non-zero expectation value of the density, pj jS (k,m), 
within one replica. The full statistics of the local static 
fluctuations that we need for the characterization of 
the macroscopic symmetries in the gel is specified by 
the order parameter field f2(k, m), encoding moments 
of arbitrary order. In particular macroscopic transla- 
tional invariance requires X^ a =o ^ Q = 0> macroscopic 
rotational invariance requires X)a=i m <* = 0> l° n g range 
orientational order for crossing angle 9 = ^ requires 
m a = IM with I e Z. 

Let us now encode these physical expectations in a 
variational Ansatz: The liquid state is characterized by 
f2 sp (k, m) = S^qS^q. It becomes unstable at a criti- 
cal cross- link density, where a percolation transition 
takes place and a macroscopic cluster is formed contain- 
ing a finite fraction of the polymers. To account for the 
fraction of localized particles Q in the percolating cluster 
coexisting with mobile particles (fraction 1 — Q) in finite 
clusters, we make the following general Ansatz for the 
expectation value of the order parameter: 

ft sp (k, (p) = {l- Q) S k6 + Q w(k 2 , <f>) 5 0)E n =o kQ (22) 

The first term describes the sol phase which is charac- 
terized by perfect spatial homogeneity and orientational 
isotropy. The second part describes the amorphous solid 
which is macroscopically translational invariant as re- 
flected in S j2 n _ k a ■ The details of the amorphous solid 
phase under consideration have to be implemented in the 
order parameter w(k , ip). 

Plugging the generic form of the order parameter (22) 
into the replica free energy (14) we get 



(23) 



(1 - Q) ^ ( cxp ]T. / ^ 6 ^ o k » A(0 - 0', 0) W (k 2 , <fi) I ^5^' - 0(a))} 



yn 



Hi 



n+l 



We expect the gel fraction to be small close to the gela- 
tion transition and thus expand the log-trace contribu- 
tion of the free energy in Q. The terms linear in Q cancel, 
as they should for the expansion around Q — to be jus- 
tified. 



A. M-fold orientationally ordered amorphous solid 

1. Hard cross-links 

Assuming that the thermal fluctuations of the poly- 
mers are to a certain degree suppressed by their stiffness 
and a sufficient number of cross-links has been formed, 
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long range oricntational order may be present as sketched 
in Fig. 2 and Fig. 3. The Ansatz for the M-fold orien- 
tationally ordered amorphous solid favors M preferred 
orientational axes separated by angles ^j. A simple way 
to incorporate this symmetry is the following Ansatz for 
the replica order parameter 



w(k,<p) 



(24) 



Here Iq denotes the modified Bcssel function of the first 
kind; it ensures the proper normalization. 

For k = 0, the above order parameter is a probability 
distribution and thus specifies the local orientational or- 
der completely In experiment on the other hand one has 
access to low order moments only. The simplest physi- 
cal order parameter being sensitive to the degree of long 
range M-fold oricntational order is given by 



S 



1 N 

-Y 

i=l 

[(cos(Atyi(s)) 



cos(M^i(s)) 



(25) 
(26) 



where in the second line, we replaced the average over 
all monomers of the system by the disorder average of 
one arbitrary monomer. This expression allows to relate 
Sm to our effective one-particle theory (. . .) denotes the 
thermal expectation value of the system for a given in- 
stance of disorder C. In order to establish the connection 



between the variational parameter r\ and the physical or- 
der parameter Sm we evaluate (26) using (24) and find 
that Sm is in leading order proportional to 77. 

We plug the Ansatz (24) in (23) and perform the k 
summations and ip integrations. If we have chosen M ac- 
cording to the symmetry considerations of the previous 
section the constraint A drops out. We are left with the 
free energy as a function of three variational parameters: 
the gel fraction Q, the spatial localization length £ and 
the degree of orientational order as measured by ;/. To 
determine these parameters we minimize the free energy. 
The resulting equation for the gel fraction is universal 
and has been derived previously [21]. There are two so- 
lutions Q = and Q ~ ^r(/i 2 — 1) which implies that the 

phase transition from sol to gel takes place at /j 2 = 1. 

In order to determine the remaining parameters ^ and 
77 we assume them to be small near the transition and do 
a Taylor expansion of the free energy where we leave out 
terms which do not depend on £ or r\ because these terms 
are irrelevant for the variation of T . More precisely, it 
will turn out that ^ oc Q, being thus small close to the 
transition, and we will keep contributions up to order 
l/£ 2 in our expansion. As for the variational parameter 
77, we will find that it actually jumps from to a finite 
value, i.e that there is a first order oricntational transi- 
tion. We expand all the same in r\ and obtain at least a 
qualitative picture. Details on the calculations are given 
in Appendix D. 

Since our Ansatz is replica-symmetric it is straightfor- 
ward to extract the part of T linear in the replica index 
n and we find 
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n 2 Q 2 



1 L " 



Ip-fC— ) 

,M 2 L 



1M,V 



J P 



~ ,M 2 L 



31 



(27) 



T ( i,, M "^s\ 7 ,/ n~,M 2 L^\ 31 a( n = M 2 L N 



Note that for M > 1 there are no terms coupling the by 
orientational order specified by 77 and the positional lo- 
calization specified by £ 2 . Such terms appear in higher 

order but will not be considered here because we restrict , \ 1 1 ■ / \ -,1 

ourselves to the vicinity of the gel point. The persistance 6 30 3 

length L„ and the localization length f are both rescaled 7 , , 8 ~ . . ,48 

, ° , , T , T h{x) ~ 1 x , h(x) ~ 1 x 

by the contour length L 01 the polymers. 21 93 

l(x) ~ x . l(x) ~ x 

K ' 6 10 ' K ' 3 5 

?/ \ 125 61 

l(x) x 28 

w 72 30 y ' 

The functions g, h, h, h and I, I, I go to zero for large ar- 
gument and are for small argument approximately given For the definitions of these functions see Appendix F. 
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Minimizing T for M > 1 with respect to £ 2 yields 
1_ J2_jA?_ = fSQ 

Hence the filaments are localized as soon as a percolat- 
ing cluster of cross-linked chains has formed and the gel 
fraction is finite. The localization length is independent 
of M and its scale is set by the radius of gyration R g of 
the filament which is ~ L 2 for stiff chains and ~ L for 
random coils. 

For finite polymer flexibility the incipient gel has 
no long range orientational order. Increasing the number 
of cross-links beyond /i 2 = 1 we find a first-order tran- 
sition that corresponds to a minimum of the free energy 
at nonzero 77, which is first metastable and eventually 
becomes the global minimum (see Fig. 4). The criti- 
cal values of /j, 2 for the first appearance of a metastable 
ordered state (//?), the first order phase transition (//?,) 
and the disappearance of the metastable disordered state 

(/i§) are shown in Fig. 5. In the limit -j- >• 00 the phase 

transition coincides with the sol-gel transition, whereas 
for more flexible filaments higher cross-link densities are 
required. For rather stiff polymers we find that the tran- 
sition takes place at 

[i\ ~ 1 + 0.33 M 2 L/L P . (30) 

In addition, Fig. 5 shows the dependence of r) c , which 
is the value of the variational parameter r\ at the phase 
transition, on the polymer flexibility. There is a ten- 
dency to a lower degree of orientational localization for 
higher values of the polymer stiffness. This behavior is 
qualitatively different from the behavior predicted for the 
lyotropic nematic ordering of partially flexible rods by 
Khokhlov and Semcnov in the corresponding range of 
flexibility [5]. However, a direct comparison of the two 
models is not feasible. In our case, orientational order- 
ing is due to the cross-links and the critical density of 
cross-links fj, 2 of the transition decreases for increasing 
polymer stiffness. There are thus two competing effects: 
stiffcr polymers should lead to a stronger orientational 
localization of the polymers whereas the smaller number 
of cross-links should have the opposite effect. It seems 
that the lower number of cross-links plays the dominant 
role in the dependence of rj c {L/ L p ). 

We point out that the cross-linking angle 9 = 2ir/M 
enters the free energy as a rescaling of the polymer flex- 
ibility L/Lp through the parameter M 2 . This implies 
that the higher the value of M (the smaller the angle), 
the higher the polymer stiffness required for the transi- 
tion at a given cross-link density. The reason for this 
scaling is that within our mean-field description we are 
dealing with one single chain in an effective medium. M- 
fold order has to be propagated along the chain from 
one cross-link to the next and the M 2 scaling reflects the 
properties of the WLC in the calculation of the corre- 
sponding correlator. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

FIG. 4: Plot of the orientational free energy of the 3-fold sym- 
metric case for L/L p = 0.05 and /j 2 = 1.1470, 1.1518, 1.1535 
(continuous, dashed, fine-dashed line, respectively). 




FIG. 5: Variational parameter of the incipient ordered state 
r/c and cross-link densities pf, /-t| and fi 2 where a metastable 
ordered minimum appears, where it becomes the global mini- 
mum (phase transition) and where the minimum at zero cor- 
responding to the disordered state becomes unstable (plot for 
M > 1). Note that always [i\ < /j 2 , < /i§ for finite stiffness. 



For a first-order transition the expansion in rj 2 is not re- 
ally justified and can only give qualitative results. Even 
worse, the expansion of T in r] 2 exhibits an oscillating 
behavior: the orders 2, 6, 10, . . . provide a stable orienta- 
tional free energy (for large 77) in a region around /i 2 = 1 
and -y- = 0, but the orders 4, 8, 12, . . . are always unsta- 
ble in the above region because they diverge asymptoti- 
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cally to minus infinity. However, considering the purely 
orientational free energy 



fo 



In f 



If / In ( f 
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(31) 
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it is obvious that for any finite k the Gaussian part is 
larger than the log-trace contribution as long as the cross- 
link density n 2 does not become too large. So, as long 
as we restrict ourselves to a region close enough to the 
sol-gel transition asymptotic stability is garantied and 
a qualitative picture can be obtained by truncating the 
expansion at order 6, f or even higher. 

For the polar case, i.e M = 1, there arc additional 
terms ins the free energy (27) which couple spatial and 
orientational part. At first sight it might be tempting to 
argue that they can be neglected close to the transition 
because they are proportional to ^ ~ Q. As it turns 
out this is not correct: the orientational transition for a 
given polymer stiffness is shifted to significantly higher 
values of ji 2 . But all the same, the qualitative picture of 
the transition is still valid. 

In order to analyze the orientational transition we first 
calculate the stationarity equation with respect to l/£ 2 

^ (32) 



'9{h)+v 2 Kh)-h i m+h i m 



and plug it again into the free energy. Keeping contri- 
butions up to order rf we obtain again a power series in 
rj 2 . The transition scenario that we found for M > 1 is 
still valid, but the transition takes place at larger values 
of fx 2 . For rather stiff polymers we find approximately 

H% ~ 1 + 0.94 L/L p (33) 

and the critical /i 2 is considerably larger than what we 
would obtain from (30) without the coupling terms. 

At first sight, it seems a bit curious that the case M = 
1 is set apart by its coupling term. This is, however, 
only due to the low order expansion and the coupling 
terms for M > 1 have not appeared yet. Because of 
the higher symmetry of the orientational contributions, 
only higher order spatial contributions may lead to non- 
vanishing coupling terms. 

If long range orientational order could exist only for 
rational values of 9/2ir, then it would be inaccessible in 
experiment. We show in the next section that the long 
range ordered states discussed above are also present if 
we introduce crosslinks that favor given crossing angles 
on average only. 



2. Soft cross-links 

Which changes do we expect when using soft cross- 
links (5) that do not rigidly fix the intersection angle to 



one particular value, but allow for fluctuations about a 
given mean direction? Soft cross-links surely will make it 
harder to establish long range order, so the first expecta- 
tion is that the orientational transition will need a higher 
cross-link density fi 2 to take place. In particular, the be- 
haviour in the limit of stiff rods will change qualitatively: 
In the case of hard cross-links the macroscopic network is 
a completely rigid object, even at cross-link densities just 
above the critical value of [i 2 — 1, whereas in the case of 
soft cross-links there are still the degrees of freedom of 
fluctuations around the preferred directions of the cross- 
links left. Instead of approaching a combined spatial and 
orientational transition right at fi 2 = 1 upon increasing 
the polymers stiffness, it appears possible that for soft 
cross-links the long range order transition will take place 
not at [i 2 = 1 but at a higher cross-link densitiy because 
the networks needs to be stabilized. However, in the limit 
of stiff rods and hard cross-links we should recover the 
combined transition right at /i 2 = 1. 

We now discuss the modifications of the free energy 
due to the soft cross-links: Keeping terms only up to 
order Q 3 as before, we find that in the case M > 1 there 
is still an effective decoupling of spatial and orientational 
transition. The spatial part does not change at all, but 
the purely orientational contribution becomes 



for. so ft 

i?o 2 f 



(34) 
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The soft cross-links give rise to additional factors ^Mgili 
. . . °^ 7 ^ 

that equal 1 in the limit of hard cross-links 7 — > 00. 

They appear linearly in the Gaussian and quadratically 
in the contribution of the log-trace term. Their range 
lying between and 1 it is clear that the log-trace contri- 
bution is smaller with respect to the Gaussian so that the 
transition occurs at a higher cross-link concentration (as 
compared to the case of hard cross-links). This is con- 
firmed by a numerical analysis of the Taylor expansion 
up to 10th order in rj. 

The effect of soft cross-links on the orientational phase 
transition is illustrated in Fig. 6 for M = 3: First of all, 
it shows that for finite 7 the phase transition takes place 
at a finite distance from /j, = 1, even in the limit of 
stiff rods. Moreover, comparing the curves correspond- 
ing to increasing values of 7, i.e. to harder and harder 
orientational cross-link constraints, the curves converge 
towards the solid curve at the bottom that was drawn for 
perfectly hard cross-links as they should. 

The case M = 1 involves additional coupling terms 
between spatial and orientational parameters that need 
to be calculated. But as before, they don't lead to a 
behavior that differs qualitatively from what we found 
for M > 1. 
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FIG. 6: Plot of the critical /i 2 (L/L p ) where the phase tran- 
sition takes place for M = 3 and stiffness parameters 7 = 
80, 160, 240, 320 and hard cross-links. 



local preferential axes by e we get 

t A n 



a=l 



(35) 



Altogether the order parameter for polar filaments in the 
glassy state then reads 



w(k, tp) 



1 



TO 



(36) 



What is the physical order parameter and how does it 
scale with the parameter 77? 

Locally, we have for each localized polymer segment 
a polar moment (tj(s)) ~ rj 7^ 0, but because of the 
orientational disorder it vanishes globally. We thus need 
an Edwards- Anderson like order parameter and choose 



1 N 

-Y 



(1*00) • (u(s)} 



(37) 



As for the M-fold order parameter in the preceding sec- 
tion, we relate the variational parameter rj to q and find 
that in lowest order 



In the preceding section we found that the higher M 
the more cross-links are needed to get into the long range 
ordered phase. This holds true for soft cross-links, too, 
as we checked numerically. The related calculations are 
presented in Appendices D and F. 



B. Statistically isotropic amorphous solid (SIAS) 

In the preceding two sections, we found that (given a 
suitable cross-linking angle 9) there is a phase boundary 
/i 2 (jj-) above which long range orientational order be- 
comes possible. For lower cross-link density, long range 
orientational order vanishes. But the positional local- 
ization of the polymer segments that takes place in the 
macroscopic cluster is always accompanied by orienta- 
tional localization. The corresponding alternative to long 
range order are glassy states, where the average orienta- 
tions of the polymer segments are frozen in random direc- 
tions, so that isotropy is restored on a macroscopic level. 
We expect such states for WLCs with a small persistence 
length, such that the order induced by a cross-link can- 
not be sustained along the contour length up to the next 
cross-link. This will be particularly severe, if M is large. 

Frozen orientations for polar filaments are described 
by the distribution (18) for a single site. However the di- 
rection of localization varies from chain to chain, so that 
averaging over the whole sample implies an average over 
the locally preferred orientation ipo assuming all direc- 
tions to be equally likely. For convenience we introduce 
the unit vectors u Q = (cos (p a , sin ip a ) and denoting the 



q = 



o(v 4 ) 



(38) 



We plug the Ansatz (36) for the replicated order pa- 
rameter into the saddle point free energy and keep only 
terms which depend on £ or rj. We obtain the simplest 
non-trivial free energy by including terms up to second 
order in i] 2 . 



n 



/' 



ln 77 772.9(^) 



6 e ' ^ 4£ 2 3 v L p 



(39) 



-^A( 7 ) ( 1 - M 2 A( 7 )M^) ) + ^^A( 7 
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rj L 



8^ 



L 



where we, have introduced the shorthand notation 
A( 7 ) := T^XtJ' mnc ti° ns 9,h,l are given by (28) 

in the M-fold section. Minimizing the above free energy 
with respect to £ 2 , yields qualitatively the same result as 
for the long range ordered state, namely 



L 2 2 n 2 Q 



e 3 5(£) 



(40) 



The orientational part shows a behavior different from 
the M-fold case: The stationarity equation with respect 
to rf gives rise to 



L 2 ^) 



(41) 



The coupling terms implies a non-zero value rj 2 , i.e. ori- 
entational localization, as soon as positional localization 
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FIG. 7: Phase diagram in the plane of cross-link concentra- 
tion /j 2 and polymer flexibility L/L p . To the right of the 
vertical dotted line, SIAS order becomes stable, wheras M- 
fold order appears to the right of the continuous tilted line 
which depends on M. 



the isotropic sol as well as below the free energy of the 
SIAS. Hence, we expect it to win as soon as it appears, 
even though we cannot do a complete stability analysis 
beyond the variational Ansatz. 

We find that the appearance of long range order is 
pushed to higher values of /_t 2 as the constraint for the 
crossing angle is softened. The same effect is observed 
for increasing flexibility of the polymers, - because it 
becomes more difficult to sustain the orientation of the 
polymers and increasing M. 

When reading the phase diagram, we should keep in 
mind that the parameters /x 2 and L/L p are not thermo- 
dynamic ones like a temperature or a chemical potential 
because changing either of them changes the disorder en- 
semble, too. This means that two points in the phase 
diagram correspond effectively to two different systems. 



sets in. Hence glassy orientational order is enslaved to 
positional localization and the orientational order param- 
eter grows continuously at the gelation transition. Vary- 
ing 7 we see that the softer the cross-links are, the smaller 
is also the variational parameter n, i.e. the degree of ori- 
entational order. 

Note that, in contrast to the M-fold ordered case, the 
limit n — > leads in the SIAS free energy to an extra 
minus sign for all the orientational contributions because 
of the coupling of the replicas in the corresponding order 
parameter (36). As a consequence, we have to maximize 
the free energy with respect to 7/ instead of minimizing 
it as it is well known from spin glasses [27]. 

C. Phase diagram 

The results of the previous sections can be summarized 
in a phase diagram, presented in Fig. 7. The control 
parameters are the cross-link density measured by /i 2 , 
the polymer flexibility measured by L/L p , and the cross- 
linking angle 8. Irrespective of the cross-linking angle, 
independent of the stiffness of the filaments and of the 
softness of the cross-links there is a continuous gelation 
transition accompanied by random local orientational 
ordering (SIAS phase) at the critical cross-link density 
/i 2 = 1. This glassy ordering has been encountered pre- 
viously for randomly linked molecules with many legs 
(p-Beine) [20, 22]. The free energy of the SIAS is above 
the free energy of the sol, which however is unstable be- 
yond /i 2 = 1 and hence, is not available in this region of 
the phase diagram. 

What is new, is the appearance of a state with long 
range orientational order, if the crossing angle 9 = 2tt 
where k, M G Z. This phase is characterized by a spon- 
taneous breaking of the rotational symmetry. For sensi- 
tive cross-links the symmetry of the orientational order 
is M-fold, for unsensitive cross-links and odd M the re- 
sulting phase has 2M-fold symmetry. The free energy of 
the long ranged ordered state is below the free energy of 



V. CONCLUSIONS - DISCUSSION 

In this paper, we studied the role of the cross-linking 
angle in the formation of orientational order in random 
networks of semiflexible polymers in two dimensions. We 
have used a variational Ansatz to map out a phase di- 
agram. Besides a statistically isotropic amorphous solid 
(SIAS) we find more exotic gels with random positional 
order coexisting with long ranged orientational order, 
whose symmetry is dictated by the crossing angle. In 
analogy to liquid crystals with long range orientational 
order and thermal centre of mass motion like in a fluid, 
these gel phases might be termed "glassy crystals" - with 
long range orientational order and frozen in random po- 
sitions like in a glass. It is interesting to note that the 
tctratic ordering, which corresponds to M = 4 in our 
system, has been predicted and/or observed in a variety 
of physical systems with quite different constituents and 
underlying physical mechanisms [14, 30-32]. 

Because of the peculiarities of two dimensions, we ex- 
pect fluctuations to affect positional localization [33, 34]. 
It is known, in the context of 2d defect-mediated melt- 
ing, that at finite temperature positional order can only 
be quasi-long ranged whereas orientational order can be 
truely long ranged [35]. A study of the corresponding 
phenomena for our positionally amorphous and orienta- 
tionally ordered system is a very interesting direction for 
further investigation. 

We point out that the finite bending rigidity is an es- 
sential ingredient of our model. It allows the effective 
decoupling of positional and orientational degrees of free- 
dom close to the gelation transition. In the case of in- 
finitely stiff polymers on the other hand, a rigid cross- 
link would automatically fix both position and orienta- 
tion. The nature of the emerging network is an interest- 
ing problem which goes beyond the scope of this paper. 

Another possible extension of our work concerns more 
elaborate variational Ansatze probing the appearance of 
combined M-fold and glassy orientational order. Here, 
the chains in the gel fraction are assumeed to be orien- 
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tationally localized in preferential directions which vary 
from chain to chain but macroscopically average in an 
M-fold pattern. 

The three-dimensional generalization of our model has 
to deal with the fact that a finite cross-linking angle be- 
tween two wormlike chains prescribes a cone and not a 
plane. If we want to describe cross-links with torsional 
rigidity, we need to go beyond the simple WLC model 
and use the helical WLC [36]. 

In this work, we focused on orientational and glassy 
order in networks of scmiflcxiblc polymers mediated 
through appropriate cross-links. We assume that the 
cxcludcd-volume interaction is such that it supports a 
macroscopically translationally and rotationally invari- 
ant liquid which, upon cross-linking, may give rise to or- 
dered networks. Although this assumption is mathemat- 
ically consistent and facilitates the analytical treatment 
of our model, it may be challenged in experimental real- 
izations where the excluded volume may lead to lyotropic 
alignment in dense systems. In a future extension of our 
model, one may envisage adding a Maier-Saupe aligning 
pscudopotential as in Refs. [6, 7] and studying its inter- 
play with the cross-link induced interaction. 
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Appendix A: Evaluating expectation values - the 
WLC propagator 

The wormlike chain propagator G(ip, s±; ip f , s 2 ) quanti- 
fies the probability that the tangential vector of monomer 
Si points into the direction ip provided that the tangen- 
tial vector of monomer s 2 points into the direction ip': 



G(ip,s 1 ;ip / ,s 2 ) 



5(p - ^( Sl )) S(p' - ^(s 2 j) 



(A.l) 



7-twLC 



Here, Af denotes the normalization and if>(s) the angle 
corresponding to the unity vector t(s). In principle, the 
path integral includes all the monomers from to L, but 
the monomers that lie not between s\ and s 2 do not affect 
the result and can be integrated out. In order to perform 
the remaining path integral we write down a discretized 
version of the above expression replacing the continuous 
degrees of freedom by a finite number I of them such 
that ti corresponds to t(s 1 ) and that t; corresponds to 
t(s 2 ). The distance e between neighbors is determined by 



le = \s2 — s±\. Expressing the integral and the derivatives 
of the wormlike chain Hamiltonian by their discretized 
versions and calling the normalization constant for the 
discretized path integral Af e we arrive at [37] 



G, 



exp 



2-1 

tE 



(1 - COs(lfii - ipi+l)) 



We want now to perform the integrations over the tpi. For 
that purpose it is convenient to decouple ipi and <Pi+\ by 
means of 



E 7 <» ei9 



q— — oo 



where the I q (a) denote modified Bcssel functions. Per- 
forming the integrations we arrive at 



1 



E 



%(-) 



2-1 



jq(<pi-tpi) 



q— — oo 
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and integrating over jp- and ^ we find that the nor- 
malization is given by Af e = (exp(— K/e)Io(K/e)) 1 . 
Last, we take the limit e — > keeping le = |si — s 2 \ 
constant. It is thus possible to express the modified 
Bessel functions by the asymptotic expansion I q (a) ~ 

^= (l — ^g" 1 + . . . ) and the propagator converges to 



G((p,s 1 ;(p',s 2 ) 



q— — oo 



e 2 - 



■<? 2 l s i- s 2| p iq(v>-v>') 



(A.2) 



Calculating a general 2-point correlation function of the 
(real valued) observables 0\ and 2 at positions s\ and 
s 2 respectively by means of the above propagator we find 



Hi 



Oi(V(*i)) 2 (1>(82))) (A.3) 

e-^« a i-i-»i d 1 (q)6;(q) , 



E 

q— — oo 



where (D\ and 2 denote the Fourier transformation of 
the observables. We learn from this expression that only 
the Fourier components to the same q couple to each 
other and that the rate of the exponential decay of cor- 
relations along the filament scales with q 2 . 



Appendix B: Mean-field replica free energy and the 
saddle point equations 

In order to obtain the disorder averaged free energy 
[F] by means of the replica method we need to calculate 
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the disorder averaged n-fold replicated partition function 

1 " 



m -En 




M\ \2Ny 



M < 
e=l I 



(B.l) 



n+l 



where vf(s) and ipf(s) denote the position vector and 
angle of orientation of segment s belonging to polymer i 
inside the crth replica. For sensitive cross-links 9 a equals 
always the single crossing-angle 0, i.e. the normalization 
y = 1, and we can omit the summation over er, but in 
the unsensitive case 9 a takes the two values 9i = 9 and 
62 = 9 + tt and so, y = 2. 

Observing that the formula factorizes in the cross-link 
index e it is possible to perform the sum over the number 
of cross-links M that leads to an exponential function 

[Z n ] cx /exp^J^J siTitf-W))* 



xA(i( s )-^( s '),9) 



H 



n+l 



where for sensitive cross-links the function A is defined 
as 



(B.2) 



and in the unsensitive case on the other hand as 

if e 7lZ a cos(i/)°-e) e 7l] Q cos('i/j Q -(e+7r)) 

A " w ' g)s 2l m + — w> — 

(B.3) 

After the disorder average, all the sites, i.e. all the 
polymer segments, are equivalent but still appear explic- 
itly, coupled by the cross-linking constraint A and the 
delta functions. Expressing the delta functions in Fourier 
space we can rewrite our formula in terms of the quantity 
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A' 



Q(k,fh) = -J2 /'' ^ ~ - 
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Using this definition and writing T~L ev explicitly the repli- 
cated partition function reads 



[Z»] ex / ex p(|^^A A |Q(k, 



exp 
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In section II the parameter A 2 was introduced in its most 
general form depending on |k| and m, but as it turns 
out a constant is sufficient for our purpose and allows in 
the following for a more compact notation. Because of 
the symmetry |Q(k, to)| 2 = \Q(— k, — to)| 2 it is only the 
real part of A„ that contributes and we thus redefine the 
sensitive kernel A s accordingly as 



A, 



•(7) 



■ cos 



(B.4) 



For unsensitive cross-links A u> m equals zero if the sum 
m a is not even, but takes otherwise the values of 
the sensitive kernel in (B.4). 

We are now going to rearrange the contributions of 
intra-polymer repulsion and Deam-Edwards distribution 
into contributions belonging to the following subsets of 
the space of (n + l)-fold replicated vectors k: 

• The 0-replica sector (ORS) consisting only of 

• The 1-replica sector (IRS) including all vectors of 
the form k = (0, . . . , k Q , . . . , 0) where k" ^ 

• The higher-replica sector (HRS) containing all the 
k where wave vectors in at least to replicas are non- 
zero, i.e. there are a 7^ (3 <G {0, 1, . . . , n} with k a ^ 
and ^ 

In the following we denote the sum over ORS and HRS 
by an d the sum over the IRS as ^ k . For the IRS 
we obtain the new kernel 



A Q 



A 2 iV 
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/3/« = l 
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Using Q Q (k, to) as shorthand notation for Q(k, to) when 
only the wave vector k Q in replica a is non-zero the ex- 
pression reads 



[Z n ] cx (exp(-7V^£ A" |Q Q (k,m)| ; 



a=0 



n+l 



In order to decouple the sites we are now going to apply 
a Hubbard-Stratonovich transformation. The symmetry 
Q(k,rh) = Q*(—'k,—m) and the analogous relation for 
Q Q (k, to) have to be reflected by their corresponding 
fields J7(k, to) and Cl a (k, to). After the transformation 
[Z n ] can be written as 

[Z n ]<x Jv{n a ,Sl}cxp (-Arj"({fr,ft})) (B.6) 

where the integrals over the complex fields are meant to 
be integrations over real and imaginary parts separately. 
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The replica free energy J- is given by 



Appendix C: Stability of the 1-replica sector 



T 



Ct=0 



EE k Ei^i 2 + ^rE k Ew 2 

:0 fh fh 

ln/exp(2i^^ k ^V / AK(0« / 



(B.7) 



| — ik Q r Q (s)| exp | — imip(s)\ 

fa 

| — ikf(s)| exp | — im-0(s)| 



exp 

2 



exp 



n+l 



Using the saddle point approximation we replace the inte- 
gral (B.6) by its maximal contribution. The correspond- 
ing values of the order parameter fields have to fulfill the 
self-consistency equations that arise from 



dm{k 0l m ) 







and 



d^sQ (k , toq) 



0. 



Note that there are many ways to perform the HS trans- 
formation, each of them leading to a different replica 
free energy J- . The resulting saddle-point replica free 
energy on the other hand is always the same as it can 
be checked easily. As a consequence, we are free to rede- 
fine the fields for later convenience by doing the replace- 
ments 5ftfi(k, to) — > \/Am 3?f2(k, to) and for the imagi- 
nary part accordingly. The advantage of this transfor- 
mation is that the saddle-point fields Ogp(k, to) are now 
directly related to an expectation value of Q(k, to) as pre- 
sented in the main part, equation (12). This connection 
between Q and its corresponding field is derived by 
means of an external field that couples to Q. Performing 
then the Hubbard-Stratonovich transformation and tak- 
ing the logarithmic derivatives with respect to real and 
imaginary part of the field on both sides of the equation, 
i.e. for both the microscopic and the field theoretic rep- 
resentation of our theory, results in the desired relation 
between Q and Q. 

The ORS/HRS part of the free energy in terms of the 
new fields reads then 



2V 



■E k E A N 



(B. 



ln (-p(£rE k E A ^ 



exp 



| — ikf (s)| exp | — irhip(s) | 




n+l 



The IRS part is treated in the next section. 



The fields of the 1-replica sector, ^ Q (k, to), can be sub- 
divided into the fields where to is such that only the en- 
try to" is non- vanishing and those corresponding to the 
other possible values of TO. In the former case the fields 
describe M-fold orientationally symmetric density fluc- 
tuations with wave vector k in replica a and have thus 
a clear physical meaning. The other fields break replica 
symmetry in the sense that they describe density fluc- 
tuations e.g. in replica a accompanied by purely orien- 
tational fluctuations in other replicas. These fields are 
unphysical and need thus to equal zero. 

In order to study the stability of the IRS with respect 
to fluctuations in the physical fields p a (k, to) it is suffi- 
cient to study stability for one replica only. The expan- 
sion of the corresponding free energy reads up to second 
order 



RS 



+2 



5>(k,m)| 2 

m 

^2 P( k l' TO l)p( k 2 ; m 2) 



(C.l) 



e i(kiri+k 2 r 2 ) e i(mii/ , l+'n2V , 2) 



C(ki,k2;mi,m2) 

with r x := r(s x ) and ifj x := ip{s x ). Choosing A 2 > p 2 the 
kernel A m = 4p — ^~ cos(toA) is positive because 
j-ffi) < 1- R is thus possible to redefine the fields by 
introducing p(k, to) := -\/A m p(k, to) without changing 
the stability. The corresponding free energy is given by 



RS 



E k E 



A 2 -^ 2 ^gcosM0 



|p(k,TO)| 2 (C2) 



+2^ p(ki,mi)p(k 2 ,m 2 )C(ki,k2;TOi,TO2) 

z 'Kl,K2 ' ' 



Let us now consider the two contributions to Tirs sep- 
arately: It is evident that the first term is a positive 
quadratic form provided that A 2 is large enough. If we 
consider stiff rods as the limiting case of very unflexi- 
ble wormlike chains the matrix of the second quadratic 
form can be diagonalized by switching back from the 
Fourier modes to to the real space variables (p. Writ- 
ing r(s) = i"o + st and integrating over ro and t we find 
for C 



i(kiri+k 2 r 2 ) 



8(<pi - (<p 2 - ip 2 ) 
sin 2 (kiti/2) 



(kiti/2)' 



(C.3) 



and know thus that the second quadratic form is positive 
semi-definitc. Altogether we have thus proved stability 
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in the case of stiff rods, but we assume that this result Considering unsensitive crosslinkcr the correspond- 

will at least hold for wormlike chains with large L/L p , ing expression reads 

too. 



Appendix D: Variational free energy: M-fold 
symmetric long ranged order 

In this Appendix, we present some details on the 
derivation of the variational free energy (27) of the M- 
fold orientationally ordered amorphous solid. We insert 
the corresponding replica order parameter Ansatz (24) 
into the general replica free energy (23) and expand it in 
powers the variational parameters i and rj. 



fo = 1+nln fl+2 V 6 Mq ,2Z {ft] J fffi cos(Mq6)) 
V ^ 7 o(l) Ml) 1 

(D.2) 

i.e. for odd M the contributions corresponding to odd 
q are projected out. For M even on the other hand, we 
obtain the same result as for the sensitive cross-links. 

Note that for hard crosslinks (7 = 00) there is an al- 
ternative expression for the /„. For sensitive cross- links 
it is given by 



1. Gaussian part 



The Gaussian part reads 



fa 



2V 7 



-£ 2 k 2 



. e r lJ22=l( cos ( M V?)+cos(Mip2 )) 



The spatial part is easily computed by replacing the sum 
over the replicated Fourier variables k by an integral and 
representing the delta function in Fourier space. Per- 
forming the resulting integrations we arrive at 



fc = 



( 1 



Jo 



n + 1 V 47r £ 2 , 
For the orientational contribution f a we find 

fo 



(D.l) 



A(ifx - <p 2 , 



<fl,<P2 

1 + n In 



e vY22=i( cos ( Ml Pi)+ cos ( M 'P2)) 

gr7(cos(M(/3i)+cos(M(p2)) 



A(9?i - tp 2 ,t 



^0 (n) 



1 , 1 Lor fgW JM 

l + nln 1 + 2 > , — 



Ilirj) I Mq {l) 
(7) 



cos(Mq9) 



For 6 = ^2tt the cosine equals 1 and vanishes. In order 
to obtain the last line we replaced the three functions de- 
pending on the t^'s by their Fourier representations. Per- 
forming then the integrations lead to Kronecker Deltas 
that cancel two of the three sums over Fourier modes. 

Altogether we find for the Gaussian contribution in 
linear order in the replica index n 



fo 



(D.3) 

/ f gi7(cos(A/ipi)+cos(A/ V2 )) \ 



1 + n In 
1 + n In 

1 + n In 



Jo(2> ? cos(^))' 



3 7 ) (cos(M( V >+e))+cos(A/ V )) 



,2r,(co S (A/ V +Ma))co S (^)) 



and in the unsensitive case we have 



fo ~ 1 + n In 



1 1 



J (2r ? cos( — )) 



+J (2 ?7 cos(^- ")) 



(D.4) 



2. Log-trace contributions 



Let us now turn to the log-trace contribution that we 
will expand up to the third order in the gel fraction Q: 



fu = (j?Qfit,i + 



fit, 1 — fu,lj 



/«,3 + Vua - 3/h,i/h, 2 ) + 0(Q 4 ). 
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The first-order term gives a trivial contribution because In the second step we use f (s) = f (0) + Jjf dr t(r) and 
of transform the original path integral T>{r(s)} into an 

integral ffi+i and a path integral X>{?/>(s)} over angular 

(D.5) 

variables. Performing the r(0) integration we get a 
Kroneckcr delta setting k to zero. The last line follows 




f _^ 2 k 2 Kronecker delta setting k to zero. The last line follows 

/ "o,E« =0 k< * c x from the fact that the functions are normalized and the 



75Dq = i cos(tp a -ip"*-0) e ri J2 a cos(mc/=) ^ ^ 

X 




tf(7) TO 



two integrations with respect to ip and ip( s ) simply lead 



The second-order contribution reads 



fit; 



y2n 



i a (k*i a +k 2 3 )/2 i(kif(«i)+k a f(«2)) 



A (ft - ft,0) A (ft - $,0) 5{<p' x - Vi(si)) % 2 - ^(s 2 )) 



ViVi Jt P2<p 2 



(D.6) 

e '7lDa=l{ cos ( M Vl)+ cos ( M y2)} \ 



Before Taylor expanding the expression in the variational 
parameters, we need to perform the summations over ki 
and k 2 . We integrate over df(0) as before and get a 
Kronecker Delta imposing ki = — k 2 . There is thus only 
one summation left. We replace this sum by an integral 
and represent the Kronecker delta in Fourier space. After 
the two Gaussian integrations we find 



using the shorthand notations f a = J S2 ds t a (s) and 
ip" := ij} a {si). This expression can be expanded up to 
the desired order and the remaining task then consists in 
calculating the corresponding correlation functions. De- 
tails are given in Appendix F. 



flt,2 



1 



e 4 « 



1 



4< 2 



1 



(D.7) 



„lE«=i (cos(Apf )+cos(A<£j)) 
e vT.l = 1 {cos(M{^+^))+cos(M(4q+A^))} 



TO 



For the calculation of the third-order term, we proceed 
the same way as before and need thus to perform three 
sums over wave vectors. 



f* 1 dr t Q (t) the result reads 



With the abbreviation f" = 
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yn ^4^2 J ^3^2 y + 1 
exp • 

I* 

d iE: = i (cos(A v °)+cos(A ¥ 3°)+cos(A v °)) v J2Z = i { cos ( «Wi +A^° ))+cos(M( i/>J +A ))+cos(Af(^f +A^J ))} 



(s^f^TT ^ (frff +f2^+frf|)-E((ff) 2 + (f 2 Q ) 2 + fi Q f2 Q ))| 

l ^ \ n + a,,3=0 a=0 ' J 



' Al/3i,A</32,A</3 3 



e' 



Appendix E: Variational free energy: SIAS 
1. Gaussian part 

In the following, we present the derivation of the SIAS 
variational free energy. Let us first turn to the Gaussian 
contribution. It is given by 



/j 2 <3 2 ^ _ f 2£2 



(E.l) 



Ig"(7?) 



where u Q = (cos , sin ip a ) . The spatial part is treated 
the same way as for the M-fold case, but the orienta- 
tional contribution is slightly more complicated because 
it does not factorize in the replica index. By means of 
the integral representation of the Bcsscl function, 

,-2tt 



we get a factorizing expression and find in linear order in 
n for the orientational part of the free energy 

fo = 1 + n f In ( f A(y>i - p 2 , 0) x (E.3) 

1 



^r;(cos(ipi+)?i)+cos(< l !;2+'52) 



= l + n/lnfl + 2E|^^cos W V 

In order to obtain the last line, we simply switch from 
the angular variables to Fourier space. The expression is 
then easily expanded up to the desired order. 



2. Log-trace contributions 

The first-order contribution of the log-trace part gives 
the same (trivial) result as for the long range ordered 
case. As for the second and third order term, we again 



have to perform the k summations first, but this calcu- 
lation is done exactly the same way as before. The corre- 
sponding results are obtained from those of the M-fold 
case by simply replacing the M-fold orientational distri- 
butions in (D.7) and (D.8) by the corresponding SIAS 
distributions. 

Details on the calculation of the individual terms of the 
second order log-trace contribution are given in Appendix 
G. 



Appendix F: Evaluation of the expectation values: 
M-fold 



In our expansion we include all the terms up to first 
order in and up to sixth order in r/. 



1. Spatial part 



Let us consider the calculation of the A- term: 



-. n n \ 

±_ £ <f^>-E<r) 2 > 

a, ,3=0 a=Q J 



S2 fS2 



dr / dr' (t(r)t(r')) 



S1,S2 J SI J Si 

S2 rS2 



dr / dr' e-^l r - T 'l 



Sl,S2 •> Si J Si 

L 



(F.l) 



The first expectation value gives a contribution for a = /3 
only. Combining the two terms and noticing that each 
replica gives the same contribution we obtain the prefac- 
tor of —to. It is interesting to note that the function g is 
closely related to the radius of gyration R g . 



R 2 g ■= \ I ((r( Sl )-r( S2 )) 2 ) (F.2) 
drdr' (t(r)t(r')) 



Si ,S2 J S2 



18 



2. Coupling terms 



value. We find 



In order to calculate the lowest order coupling terms 
we expand the spatial part in first order. For 9 = jj^n 
we find 



4* 



/ 7i n 



a//3=0 " ' ~ a=0 

e T 5Za = i ( cos A Vi +™s A^J ) 



X 



lAip 1 ,A'tp 2 Iq re (7) 

c '?E" = i{cos(A/(V7 + A V 7))+co S (A/(^+Avp2))} 



We can omit the second term in the first line of the above 
equation because it is already proportional to n and will 
in the end give rise to contributions proportional to n 2 . 
Writing for the moment only the terms that are directly 
involved in the thermal expectation value, we have 



n ; 

^ / | cos ip° cos ipP, + sin ip° sin ip^, j 



xe 



'?E^ = i{co S (A/(07+A ¥ >7))+cos(A/(0j+A v J))} 



-nS M A (cos^e^+^S <Jcos^<e" (c2+ci) J> 



(*) 



+ (sinV4c" (c 5 +c ^\ (smipl,e r '^ +c ^ 



_J_ f e7 co8Av, le7 cosA V2 x (F.5) 

7o (7J JAip u Aif2 

X (cOs{%l) T ) c ^{ cos (^i+ A< />i)+ cos W>2+Aip2)} ^ 

, / u / \ ^ g (7)^+l(7) ^ 

2^ ^(iJVhi) — 727— 



I| e ^^ (<?2|si ~ T|+(<?+1)2|;i2_T|) 

+c -27-((9+ 1 ) 2 | ;s i- T l+'J 2 |^-r|)| ^ 



(sin(V>r) e" {cosW ' 1+Atpl)+cos( ^ +A¥,2)} ^ = (F.6) 



and 



e 7cosAifig7cosA(fJ2 x 
Ai£i,Aip2 

1 77{cos(i/>(si)+A(pi)+cos(j/>(s 2 )+Ai ) 02)} 



4 2 (7) 



(F.7) 



= i?(,)+2£i?(^)^e-A«"l- 



From (F.4) we obtain the coefficients of the coupling 

2 4 6 

terms proportional to ^ , |r and ^ by sampling all the 

ways to collect the corresponding power in 77 from the 

2 

different factors. For e.g. the only way to collect an 

factor r] 2 is to expand the two averages in (*) in first 
order and leave out (**). The result is 



x[[(e'K^) 

7=3 



(F.4) 



where we used the abbreviations cf := cos(-0 a (si) + 
Aipf). This expectation value is only non- vanishing for 
M = 1 because of the relation (A. 3): the Fourier trans- 
formed of cos(i/v) has contributions for the modes q = 
±1, but the Fourier transformed of exp(?7 cos(M-0(si)) 
only for q = ±ZM. The prefactor of — n ~ n(n — 1) 
stems from the number of ways to combine the indices 
a and /3 from the spatial contribution with the 7's from 
the orientational contribution. 



I 2 
2 <? 



t2t T 4 

-JfHa)-^ Ha) 



4 \ 48 7 2 



12 II 



h(a) 



16 7 n 4 



, s . 1 7i7 2 7 3 , / \ 

'5(a) + — ^ < 4 (a) 



1 r 2 r 2 r 4 r H 



(F. 



(a) 



We need thus to compute three types of expectation 



The appearing correlation functions arc defined as fol- 
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lows: 



index wc obtain 



l(x) := 

h{x) := 

h(x) := 

l 4 (x) := 

l 5 (x) := 

l 6 (x) := 

h{x) := 



7 - 80-* + c- 2x -6x + 2x 2 



2x A 



-11 - 9c 



-2x 



2e 



-3x 



18c" 



6a; 



Ox 4 V 



600a; 



9a; 4 

291 - lOOc" 2 * + 400c" 



-16e 



-5a; 



25e" 



180a: 



1 



- 653 - 1296e 



- 5:1' 



81000a; 4 

+2025c- 4a; + 324c" lte - A00e~ 9x 
1 



(F.9) 
(F.10) 
(F.ll) 

(F.12) 



7200a; 4 



297 + 400e~ 



128c 



— bx 



1260a;J 
(F.13) 



1 



1800a; 4 



iOa; 4 V 



1350a; 



+25e~ ax + 360a; 
37 + 100c~ 6x - 288c~ 5k (F.14) 



-225e _4a; + 60a; 



114-100e" 3a; + 225e" 2:z; (F.15) 



+25c" 6:r - 36e~ 5:E + 120a; J 



In the case of hard cross-links it is convenient to intro- 
duce the following functions: 



l{x) :=l&(l(x)+Al 2 (x)-h{x) 



(F.16) 



and 



l(x) := ^l( x ) + U8l 2 (x)-^-h(x) (F.17) 
+ \h{x) + l 5 {x) - 4/ 6 (.t) - 16Z 7 (z)) 



In 



3 7(cos Aipi+cos Aip 2 ) 



Aip±,Aip 2 



(F.19) 



The A(/3-intcgrations lead to a Fourier transformation of 
the soft cross-link contributions and noticing that the 
resulting expression is symmetric in q we find 



(F.20) 



Expanding the orientational free energy for hard cross- 
links in r\ and performing the s\- and S2-integrations it 
is convenient to introduce the function h(x) 



and to build of h(x) the functions 



(F.21) 



and 



h(x) = h(x) + h(2x) - -h(4x) 



h(x) = — h(x) + 3h(2x) + h(3x) 
16 

-l h (4x)-h(5x) + ±h(9x) 



(F.22) 



(F.23) 



3. Orientational contributions 

As for the calculation of the purely orientational part 
we notice first that it is factorizing in the replica index 
a. It is thus possible to take the replica limit directly by 
means of the expansion (. . .)" = 1 + rein (...)+ 0(n 2 ) 
and wc have 



for 



l + Ti / In I 

Sl,S2 



3 7(C0S Al^l+COS A 1^2) 



X 



(7) 

c j ) (cos(A/(^i+A V i))+cos(J\/(V'2+A V2 )) 



(F.18) 



The thermal average can be performed using formula 
(A. 3). Keeping only the contribution linear in the replica 



Appendix G: Evaluation of the expectation values: 
SIAS 



The calculation of the expectation values in the sec- 
ond order contributions of the SIAS log-trace part can 
be done along the same lines as above using again the 
integral representation of the Bessel function (E.2). The 
calculation of the lowest order spatial contributions is the 
same as for the M-fold case. For the purely orienta- 
tional we obtain 



flt2,or 



1) 



xe 2i 



^-"IcosfaOPi-tfa))) 
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The lowest order coupling term is given by As in (F.3) the second term in the first line of the ex- 

pression is already proportional to n and will thus finally 
1 // 1 " n n \ lead to a contribution of 0(n 2 ) that is not of interest for 

4fMl n + 1 ^°^~ n I 1 ^~"^°) 2 1 X (^'^) us - After integrating with respect to the ?9-variables and 

\V a^/3=o a=o J then with respect to the (^-variables we obtain in the end 



e 



7 YLa=i ( COS Al^"+COS A(y?2 ) 



^Ja^A^ J n W _ n jl^4 |(a) _ (G3) 

8£ 2 ^o(7) 



7,5=1 The results for the SIAS can be expressed in terms of the 

\ functions g, h and I that we already introduced in the 

x ( cos^i — #1 + Ay?f) + cos(?/>2 — $2 + A-ifz) J \ preceding section. 
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